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ABSTRACT 

The galaxy power spectrum is now a well-known tool of precision cosmology. In addi- 
tion to the overall shape, baryon oscillations and the small-scale suppression of power 
by massive neutrinos capture complimentary information on cosmological parame- 
ters when compared to the angular power spectrum of cosmic microwave background 
anisotropics. We study both the real space and redshift space galaxy power spectra 
in the context of non-linear effects and model them based on the halo approach to 
large scale structure clustering. We consider potential systematic in the cosmological 
parameter determination when non-linear effects are ignored and the galaxy power 
spectrum is described with the linear power spectrum scaled by a constant bias fac- 
tor. We suggest that significant improvements can be made when non-linear effects 
are taken into account as a power-law contribution with two additional parameters to 
be determined from the data. In addition to cosmological parameters through galaxy 
clustering, such an approach allow a determination of useful information related to 
astrophysics on how galaxies occupy dark matter halos. 

Key words: cosmology: theory — dark matter — large-scale structure of Universe 
— galaxies: distances and redshifts — Galaxy clustering: galaxies 



1 INTRODUCTION 

Through a statistical measurement of information related to 
inhomogeneities in the universe, the galaxy power spectrum 
provides a strong tool in the present era of precision cos- 
mology. In addition to information related to the primordial 
power spectrum, characterized by a sl ope and a norma liza- 
tion, through a turnover at fceq = yj2Q,mHQ{l + 2cq), the 
overall shape of the galaxy power spectrum captures cos- 
mological information related to the horizon at the matter- 
radiation equality. In addition to the overall shape, addi- 
tional features in the galaxy power spectrum allow improved 
measurements of several cosmological parameters. Similar 
to the oscillatory features in the CMB anisotropy power 
spectrum, the matter power spectrum is expected to ex- 
hibit the presence of baryons through oscillations. These 
features are associated with the sound horizon at the end 
of the Compton-drag epoch and capture information related 
to Q.mh'^ and Q.bh? (Eisenstein & Hu 1998). 

Unlike well-known oscillatory features in the angular 
power spectrum of cosmic microwave background (CMB) 
anisotropics, baryon related oscillations in the matter power 
spectrum of the large-scale structure are highly suppressed 
due to the low baryon-to-dark matter ratio. For currently fa- 



vored ACDM cosmologies, baryon oscillations in the matter 
power spectrum have amplitude variations at the level of ~ 
8% with effective widths in Fourier space of order Afc ~ 0.03 
h Mpc~^. While the ideal way to detect baryon oscillations 
is to measure the three-dimensional dark matter power spec- 
trum directly, unfortunately, there are no useful probes of 
this quantity. The effects such as the weak gravitational lens- 
ing, which trace matter fluctuations, only provide informa- 
tion related to the projected power spectrum over a broad 
window function in the redshift space (e.g., Bartlemann & 
Schneider 2001). This leads to an averaging of small fea- 
tures such as baryon oscillations such that they remain un- 
detectable even considering most favorable scenarios. 

Galaxy redshift surveys, on the other hand, allow a 
measurement of the three dimensional power spectrum of 
the galaxy distribution. Since galaxies are expected to trace 
matter fluctuations, at least in large and linear scales, the 
galaxy power spectrum has been pursued for an obser- 
vational detection of baryon oscillations. When detected, 
baryon oscillations are expected to allow significant improve- 
ments in parameter determination by breaking various de- 
generacies associated with the interpretation of CMB data 
(Eisenstein, Hu & Tegmark 1998; Blake & Glazebrook 2003; 
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Eisenstein 2002; Cooray 2002a). While expectations for 
precision parameter measurements are generally high, these 
have been tested to some extent with measured power spec- 
tra from redshift surveys such as the 2dF Galaxy Redshift 
Survey (2dFGRS; Colloss ot al. 2001) when combined with 
recent results from CMB anisotropy experiments (e.g., Per- 
cival et al. 2002). The same surveys have allowed first at- 
tempts to detect baryon oscillations, though, there is still 
no significant evidence for them (e.g., Percival et al. 2001; 
also, see. Miller Nichol & Batuski 2001). This can be un- 
derstood by the fact that three-dimensional redshift sur- 
veys, so far, lack the required volume, with sizes of order 
L ~ 2Tv/Ak ~ 300 h~^ Mpc in each dimension, to reliably 
resolve them. 

In addition to galaxy redshift surveys, a number of addi- 
tional observational efforts are either under way or planned 
to image the large-scale structure out to a redshift of a few. 

These wide-field surveys typically cover tens to thousands 
of square degrees on the sky and include weak gravitational 
lensing shear observations with instruments such as the Su- 
pernova Acceleration Probe (SNAP; Masscy et al. 2003) 
and the Large Aperture Synoptic Survey Telescope (LSST; 
Tyson, Wittmann & Angel 2000), aud observations of the 
Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1980) with 
dedicated small angular scale CMB telescopes such as the 
South Pole Telescope (SPT; Stark ct al. 1998). These sur- 
veys are expected to produce catalogs of dark matter halos, 
which in the case of lensing and SZ surveys are expected to 
be essentially mass selected (Wittman et al. 2001; Holder 
et al. 2000). While the halo number count as a function of 
redshift is a well-known cosmological test (Haiman, Mohr 
& Holder 2001), one can also consider the additional in- 
formation supplied by the clustering of halos. In particular, 
similar to the approach in galaxy redshift surveys, one can 
measure the power spectrum of clustering associated with 
these halos, as a whole, and use that power spectrum for 
cosmological studies. In addition to the shape, cosmological 
information comes from the redshift evolution of rulers that 
can be calibrated through CMB data (Cooray et al. 2001; 
Blake & Glazebrook 2003). 

The cosmological studies based on the galaxy power 
spectrum and related halo clustering information, unfortu- 
nately, is affected by non-linearities. While sources are ex- 
pected to trace inhomogcneities in mass, due to the non- 
linear evolution of gravitational perturbations at late times, 
the galaxy power spectrum at small scales departs signifi- 
cantly from the linear description corrected for source bias. 
The enhancement of power due to non-linearities erases os- 
cillatory features (Meiksin, White & Peacock 1996), while 
effects such as redshift-space distortions further contribute 
to their disappearance. Due to the cosmological significance 
of baryon oscillations, suggestions have already been made 
to detect them at redshifts of a few since the non-linear scale 
at high redshift is expected to move to smaller scales than 
today (Eisenstein 2002; Blake & Glazebrook 2003). 

In the case of current measurements, the cosmological 
interpretation is restricted only for large scales where clus- 
tering is expected to be linear. The standard approach to de- 
scribe the galaxy power spectrum involves the linear power 



spectrum scaled by a constant bias factor (e.g., Percival et 
al. 2001). While slight modifications to the bias description 
have been considered (Elgaroy & Lahav 2003), there is still 
no guidance as to how non-linearities affect cosmological pa- 
rameter measurements and what approaches can be made 
to improve such estimates. Additionally, by restricting cos- 
mological interpretation to specific range of wave numbers, 
it is likely that we have only extracted limited information 
from current data. For example, in the case of the 2dFGRS, 
analyzes have only been considered out to a wave immber 
of 0.15 h Mpc"'^, while measurements span to considerably 
smaller scales. It is clear that further studies are needed 
on galaxy clustering well in to the non-linear regime and 
to understand what information can be extracted from the 
non-Unear regime of clustering, or at least, from transition 
regime where non-linearities become important. 

The purpose of this paper is to understand the onset of 
non-linearities in the galaxy power spectrum. While detailed 
numerical simulations are useful for this purpose, here, we 
make use of the halo approach to galaxy clustering (Seljak 
2000; Peacock & Smith 2000; Scoccimarro ct al. 2001; Sheth 
& Diaferio 2001; Berlind & Weinberg 2002). The halo ap- 
proach provides a useful tool to understand certain aspects 
of clustering and is physically consistent model without any 
prior assumptions on the biasing nature of galaxies with 
respect to the dark matter distribution (for a review, see, 
Cooray & Sheth 2002). We study how cosmological param- 
eter measurements can be improved by extracting informa- 
tion from mildly non-linear scales that are currently ignored. 
In particular, we pay attention to systematics that are in- 
troduced to cosmological parameter mcEisurements when the 
galaxy power spectrum is modeled simply with a scaled ver- 
sion of the linear power spectrum. We show that parameter 
estimates are biased away from the true values and that 
these departures can be significantly reduced when the non- 
linear contribution, at least near the regime of transition 
from linear to non-linear clustering, is modeled as a power- 
law. Since parameters related to this description is deter- 
mined from data simultaneously with cosmological parame- 
ter estimates, this approach leads to a slight degradation of 
errors associated with cosmological parameters alone. The 
main advantage is that these estimates, however, remain ei- 
ther unbiased or biased at a level that is negligible. The 
suggested approach can be easily implemented with current 
studies and when interpreting results related to completed 
2dFGRS and SDSS surveys (York et al. 2000), among oth- 
ers. The suggested approach also avoids an arbitrary dis- 
tinction as to what scales correspond to the linear regime 
of clustering while also allowing the maximum extraction of 
cosmological information from data. 

The paper is organized as following. In the next Section, 
we outline main considerations related to the halo approach 
to galaxy clustering. We refer the reader to review by Cooray 
& Sheth (2003) for full details of this approach. We make 
use of the Fisher matrix approach to quantify biases and 
expected errors on parameter measurements in Section 3. In 
Section 4, we conclude with a summary. 
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Figure 1. The galaxy power spectrum in real space (solid lines) and a comparison to the clustering in dark matter (dotted lines). We 
show the 1- and 2-halo contributions to the galaxy power spectrum and the 1-halo contribution to the dark matter power spectrum. For 
reference, the dot-dashed line shows the linear power spectrum. The measured data points represent the power spectrum of the IRAS 
Point Source Catalog Redshift Survey (PCSz; Saunders et al. 2000) as measured by Hamilton & Tegmark (2002). In the bottom panel, 
we consider the transition regime between linear to non-linear clustering. As we discuss later, the transition to non-linear clustering can 
be described via a combination of the linear power spectrum scaled with a constant, scale-free, bias and an additional contribution which 
is a power-law. 
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2 HALO APPROACH TO CLUSTERING 
2.1 Real Space Clustering 

The halo approach provides a physically motivated scheme 
to understand non-linear clustering of dark matter and var- 
ious tracer populations of the large scale structure (see 
Cooray & Sheth 2002 for more details). Under this descrip- 
tion, the galaxy spectrum takes the form of 
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denotes the normalized Fourier transform of the galaxy den- 
sity distribution within a halo of mass m. The standard as- 
sumption here is that galaxies trace dark matter within halos 
such that one can utilize the dark matter distribution given 
by analytic forms such as the NFW (Navarro et al. 1996) 
profile. An improved approximation will be to use the den- 
sity distribution defined by sub-halos. Since relevant profiles 
related to substructure is still not well studied numerically, 
we make use of the dark matter density profile to describe 
the galaxy distribution within halos. In equations and |5| 
n{m) is the halo mass function (e.g., Press & Schechter 1974) 
and b\(m) is the halo bias relative to linear clustering (Mo 
& White 1996). 

An important ingredient in the halo description of 
galaxy clustering is the halo occupation number (Kauffmann 
etal. 1999; Benson et al. 2000; Berlind et al. 2003). Follow- 
ing Sheth & Diaferio (2001), for the illustration of non-linear 
effects related to galaxy clustering, we take a description 
motivated by semi-analytic models: 



iVo 

iVo {m/M^Y 



10" MQh-^ <m< Mb 
m> Mb 



where we set Mb = 4 x lO^'^ Mq/H and take a and A^o as free 
parameters. In our fiducial description, these two parameters 
take numerical values of 0.8 and 0.7, respectively. 

In addition to the mean number of galaxies as a function 
of halo mass, for the one-halo term of the power spectrum 
(equation^, we need information related to the second mo- 
ment of the halo occupation number, (A'gai(A'gai — l)\m). If 
the probability distribution function for galaxy occupation is 
Poisson, then {7Vgai(A'gai — l)|m) — {A'^gai|m)^, though depar- 
tures from the Poisson description is clearly visible in numer- 
ical data (Berlind et al. 2003) . Following Scoccimarro et al. 
(2001), we make use of the binomial distribution, matched 



to numerical data, to obtain a convenient approximation. 
The second moment is then 

{iVgai(iVgai - 1))'''' = P(m) {iVgaljm) , (5) 

where /3(m) = log yJm/lQ"h-^MQ for m < lO"ft"^M0 
and P{m) — 1 thereafter. Note that in equation the sim- 
plest approach is to set p = 2 when calculating P]^{k). In 
halos which contain only a single galaxy, however, it is nat- 
ural to assume that this galaxy sits at the center. Thus, an 
improved model for the galaxy distribution is to take p = 2 
when (A'gai(Afgai — 1)) is greater than unity and p = 1 oth- 
erwise. 

In Fig. 1, we show the galaxy power spectrum calcu- 
lated following the halo model and a comparison to the non- 
linear dark matter power spectrum, again calculated under 
the halo model; to describe dark matter clustering, one re- 
moves information related to the halo occupation number 
and the number density of galaxies from equation Note 
that in Fig. 1 and throughout, we plot the logarithmic power 
spectrum, A^(fc) = k'-'P{k)/2n^ instead of P{k) itself. We 
also show a measurement of the real space power spectrum 
from the IRAS PSCZ survey (Saunders et al. 1990) as de- 
termined by Hamilton & Tegmark (2002). In the bottom 
plot of Fig. 1, we show the transition regime between linear 
and non-linear clustering. At large scales, where the two- 
halo term dominates, the galaxy power spectrum traces the 
linear clustering with a constant bias. One can understand 
this by noting that at large scales, Ugai(fc|m) 1 and the 
galaxy power spectrum simplifies to 



J'gal(fc) « &gal J'""(fc), 



with 



fegai = / dmn{m) bi{m) 



n-gal 



(6) 



(7) 



denoting the mean bias factor of the galaxy population. 

This is consistent with the standard assumption in the 
literature that the galaxy bias, at large scales, is constant. 
The drawback with this description is that the physical scale 
below which this assumption breaks down is ill defined. It 
is generally assumed that the on set of non-linearity is at 
wave numbers ~ 0.1 h Mpc"'^, though this is only based on 
numerical simulations of the non-linear dark matter power 
spectrum. While such a scale is valid for the dark matter 
distribution, when considering galaxies, there is no reason 
why the same situation should hold. In the case of the halo 
approach to galaxy clustering, the onset of non-linearities, 
by which we assume the dominance of the 1-halo term, is 
model dependent. To illustrate this, in Fig. 2, we model the 
galaxy power spectrum as a function of a and A'o. Note that 
the transition is strongly sensitive to a such that when a is 
low, the galaxy power spectrum traces the linear clustering 
down to relatively smaller scales, ~ 0.5 Mpc"'^, while when a 
is high (> 0.8), the onset of non-linearity is near 0.1 Mpc~^. 
If (A'gailm) simply scales as No{m/Mh)" , then the galaxy 
power spectrum is independent of No, as one renormalizes 
by the galaxy number density which is also an integral of this 
quantity weighted by the mass function. In Fig. 2, we see 
slight dependences on A^o since the halo occupation number 
was modeled with two parts, one which is a constant and 
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Table 1. Marginalized Errors 







k < 0.1 Mpc-i 




k < 0.2 Mpc-i 




k < 0.5 Mpc-i 


Parameter 


Error 


Bias 


|Bias|/Error 


Error 


Bias 


1 Bias 1 /Error 


Error 


Bias 


1 Bias 1 /Error 




1.35 


0.15 


0.11 


0.22 


0.15 


0.73 


0.20 


0.29 


1.44 


ln(n„/i2) 


1.29 


0.06 


0.05 


0.31 


0.05 


0.17 


0.29 


0.08 


1.75 


mi,{eV) oc fliyh^ 


5.64 


-1.01 


0.18 


1.59 


-1.90 


1.19 


1.35 


-7.54 


5.56 


Qa 


oo 






oo 






oo 






w 


oo 






oo 






oo 






"s(%d) 


0.58 


-0.04 


0.06 


0.09 


-0.22 


2.34 


0.08 


-0.66 


8.74 


lnP*(fefi<j) = lnyl| 


oo 






oo 






oo 






In 6 


oo 






oo 






oo 











0.03 


0.00 


0.01 


0.02 


0.02 


0.76 


0.02 


0.23 


10.56 




0.06 


0.00 


0.03 


0.05 


0.06 


1.25 


0.04 


0.39 


8.38 


mv{eV) oc 




0.38 


-0.06 


0.09 


0.23 


-0.94 


4.23 


0.21 


-3.48 


16.77 


Qa 




0.07 


0.00 


0.00 


0.06 


0.02 


0.33 


0.06 


0.16 


2.92 


w 




0.29 


-0.00 


0.03 


0.27 


-0.3 


1.11 


0.26 


-2.04 


7.85 






0.018 


-0.00 


0.04 


0.015 


-0.03 


2.02 


0.01 


-0.24 


18.67 


lnP*(fcfid) 




0.021 


0.00 


0.00 


0.017 


0.02 


1.07 


0.016 


0.19 


12.00 


In 6 




0.040 


0.00 


0.09 


0.034 


0.06 


1.93 


0.033 


0.06 


1.77 


ln(nsfc2) 




0.005 


0.00 


0.03 


0.005 


0.01 


1.17 


0.005 


0.04 


6.88 






0.015 


0.01 


0.06 


0.013 


0.04 


2.84 


0.012 


0.20 


16.21 


mi,{eV) oc 




0.17 


-0.02 


0.10 


0.12 


-0.72 


5.78 


0.09 


-3.47 


38.43 


Oa 




0.05 


0.00 


0.00 


0.049 


0.00 


0.03 


0.05 


0.022 


0.44 


w 




0.15 


-0.00 


0.03 


0.17 


-0.15 


0.89 


0.17 


-0.68 


4.01 


ns{fcfid) 




0.004 


-0.00 


0.04 


0.004 


-0.01 


2.43 


0.003 


-0.05 


16.95 


lnP<i,(fcfld) 




0.005 


0.00 


0.01 


0.005 


0.01 


1.03 


0.004 


0.04 


8.92 


Inb 




0.017 


0.00 


0.06 


0.016 


0.05 


2.92 


0.013 


0.36 


26.91 



NOTES. — Parameter error calculation assume following fiducial model: Qm = 0.3, = 0.04, SIa = 0-7, h = 0.72, w = —1 and Qi^ = 0. 
We fix other interesting parameters related to CMB, r = 0.1 and T/S = 0, and assume a flat-cosmology with Qk = 0. We quote la- 
errors. The top, central, and bottom divisions represent the parameter errors from galaxy clustering, galaxy clustering with the full 3 
year WMAP data, and galaxy clustering with Planck data. A parameter error of infinity implies that the parameter affects observables 
but is not constrained due to degeneracies; in particular, the normalization and bias are all degenerate and is only broken through the 

addition of CMB data. 



another involving a power-law. In this case, the 1-halo term 
is slightly dependent on the numerical value of A'o. 



2.2 Prom Real-Space to Redshift-Space 

So far, wo have only considered the real space power spec- 
trum of galEixies. Observationally, for example in the case 
of the 2dF survey, the galaxy power spectrum is measured 
from redshift surveys and additional corrections must be 
taken into account when modeling clustering in the redshift 
space. Following Kaiser (1987), wc can write the redshift 
space fluctuation, Sg , of galaxy density field, relative to real 
space fluctuation, 5g, as 

<5^(k) = 5g(k)-^5,/i^ (8) 

where 6v is the divergence of the velocity fleld and = f -k is 
the lino of sight angle. At linear scales, one can simplify the 
relation by noting that Sg(k) = bgS(k) and 5v = f{flm)S(k) 
with /(fim) = d log (5/d logo w to obtain 

S^g{k) = Sg{k)[l + l3n'], (9) 



where /3 = f{Qm)/bg. The latter is generally measured 
from redshift surveys since it allows a determination of Qm 
through constraints on the bias factor (Strauss & Willick 
1995; Outram ot al. 2001; Peacock et al. 2001). At 
linear scales, the distortions increase power by a factor 
(1 + 2/3/3+1/5/3^), which when bg = lis 1.41 for Qm = 0.35. 
At non-linear scales, virial velocities within halos modify 
clustering properties. With a Gaussian description for one 
dimensional virial motions, we write 

5^(k) = 5ge-('="''>'/'. (10) 

Putting all terms together, under the halo approach, we can 
write the power spectrum in redshift space as (Seljak 2001) 

Pz-gai {k) = P^-g^i (k) + P^h^^i (k) where 

-Pz— gal(^) = 

/ dmn{m) <^«^'(^f^' ~ R^{ka)\ug^iik\m)\'' , 

J "gal 

P'-\ai(fc) = + + P"°(A;) , 

(11) 
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Table 2. Marginalized Errors With non-Linearities Accounted by a Power-Law 



k < 0.1 Mpc-i k < 0.2 Mpc-i k < 0.5 Mpc~i 



Parameter 


Error 


Bias 


Bias 1 /Error 


Error 


Bias 


1 Bias /Error 


Error 


Bias 


1 Bias 1 /Error 


ln(S7s/i2) 


7.33 


0.42 


0.05 


0.42 


0.03 


0.07 


0.20 


0.005 


0.02 


ln(nm/i2) 


4.88 


0.24 


0.05 


0.38 


0.03 


0.08 


0.33 


0.02 


0.06 


{eV) oc f2^/i^ 


7.27 


-0.14 


0.02 


1.74 


-0.10 


0.06 


1.50 


-0.04 


0.03 


f^A 


oo 


- 


- 


oo 


- 


- 


oo 


- 


- 


w 


oo 






oo 






oo 






' "iD V'ria / 


0.78 


-0.02 


0.03 


0.12 


-0.01 


0.08 


0.11 


-0.01 


0.09 


lnP$(fefld) = In A| 


oo 






2.36 


-0.11 


0.06 


0.94 


0.09 


0.09 


In 6 


oo 


- 


- 


1.64 


0.11 


0.06 


0.82 


0.07 


0.08 




oo 


- 


- 


3.95 


0.31 


0.08 


0.41 


0.06 


0.14 


"nl 


oo 


- 


- 


3.58 


-0.14 


0.04 


0.45 


-0.05 


0.11 




0.25 


0.00 


0.02 


0.024 


0.00 


0.05 


0.02 


0.00 


0.05 




0.06 


0.00 


0.04 


0.05 


0.00 


0.02 


0.05 


0.00 


0.04 


'iTiv (eV') oc Q,vh? 


0.41 


-0.03 


0.07 


0.25 


-0.01 


0.04 


0.23 


-0.01 


0.04 


Qa 


0.06 


0.00 


0.01 


0.06 


0.00 


0.02 


0.06 


0.00 


0.02 


w 


0.30 


0.01 


0.04 


0.28 


0.02 


0.06 


0.27 


0.00 


0.02 


ns{%d) 


0.019 


0.00 


0.04 


0.016 


-0.00 


0.07 


0.01 


-0.00 


0.01 


lnP^{kfid) = lnA| 


0.020 


0.00 


0.03 


0.018 


0.00 


0.08 


0.016 


0.00 


0.05 


Inb 


0.13 


0.00 


0.06 


0.036 


0.00 


0.03 


0.034 


0.00 


0.06 




oo 






2.27 


0.21 


0.09 


0.22 


0.04 


0.18 


"nl 


oo 






1.98 


-0.12 


0.06 


0.23 


0.03 


0.13 


ln(nfl/i.2) 


0.005 


0.00 


0.00 


0.006 


0.00 


0.01 


0.005 


0.00 


0.04 


ln(n,„h2) 


0.015 


0.00 


0.02 


0.014 


0.00 


0.02 


0.014 


0.00 


0.07 


nil, (eV) oc Qijh^ 


0.18 


-0.00 


0.02 


0.14 


-0.00 


0.01 


0.10 


-0.00 


0.01 


Ha 


0.05 


0.00 


0.00 


0.05 


0.00 


0.00 


0.05 


0.00 


0.00 


w 


0.17 


0.00 


0.01 


0.17 


0.00 


0.01 


0.17 


0.01 


0.05 


ns(fcfld) 


0.004 


-0.00 


0.01 


0.004 


-0.00 


0.01 


0.003 


-0.00 


0.01 


lnP.j(fcfld) = lnA| 


0.005 


-0.00 


0.01 
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0.11 
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oo 






1.51 


-0.11 


0.07 


0.19 


0.02 


0.10 



NOTES. — Same as Table 1. We add two extra parameters: the normalization of the non-linear contribution to galaxy power spectrum, 
A^l and its power-law slope, a^b with scale; In the case of our fiducial model and a logarithmic power spectrum, the two parameters 
have numerical values of 65.7 and 3.04, respectively. Note the reduction in parameter biases compared to Table 1. 

where the virial radius is related to the halo mass through 
Mvir — 47rApr2ir/3 and 7 is a free parameter that we 
use here to demonstrate how uncertainties related to virial 
motion affect the predicted non-linear contribution to the 
redshift-space power spectrum of galaxies. 

Even though peculiar velocities increase power at large 
scales, virial motions within halos lead to a suppression of 
power. In figure 3, we show the power spectra of dark mat- 
ter and galaxies in redshift space with a comparison to to 
that of the real space. While clustering power is increased 
at large scales, both for dark matter and galaxies, the power 
is reduced substantially at scales where the 1-haIo term of 
the power spectrum dominates. This is easily understand- 
able since the virial motions within halos lead to decrease in 
a power when averaged over all angles. 

As a comparison to observable data so far, to model 
the 2dFGRS redshift space power spectrum, we move from 
theory space description to the observable space by convolv- 



with 



and 



dmn{m) bi(m) ^'^^"'^"^^ i?i (fco-)ugai(fc|m) , 

Wgal 



fi^m) / dmn{m) bi{m) Ri{ka)u{k\m) , (12) 



Rp{a = kay/p/2) = 



/tt erf (a) 



(13) 



2 a ' 

for p = 1, 2. In equation lllH . figai denotes the mean number 
density of galaxies (equation I^J. 

In order to model velocities within virialized halos, we 
make use of a description based on isothermal spheres and 
in reasonable agreement with measurements of virial veloci- 
ties within halos in numerical simulations (Sheth & Diaferio 
2001). If Gvir denotes the one dimensional velocity disper- 
sion of particles within a halo, then 



: ^ 



Gm 



(14) 
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Figure 2. The dependence on the real space galaxy power spectrum on model parameters related to the halo occupation number. In 
the top panel, we show variations related to a (with A^o = 0.7) and the bottom panel shows variations related to A^o (with a = 0.8). 
As shown, the slope of the mean number of galaxy as a function of halo mass strongly determines the scale at which 1-halo term, the 
non-linear contribution, dominates the total power spectrum. 



ing the theory power spectrum over the 2dF survey window 
function: 



where a = 0.00342 and h = 0.00983 (Elgaroy, Gramann & 
Lahav 2002). The observed power spectrum is then given 
by 



^2dF W = 1 ^ ^ (fc/fc)4 (15) P^aAk) = J d^k'P._gai(|k - k'lWiMk') , (16) 
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where W is the normahzed window function such that 
/ d^k Wi^y{k) = 1. 

In Fig. 4, we show predicted redshift space galaxy power 
spectra convolved with the 2dF window function. As a com- 
parison, we also plot the 2dFGRS power spectrum from Per- 
cival et al. (2001). Note that we have not attempted to model 
the 2dF data by varying parameters of either the cosmologi- 
cal model or the halo description, but rather we show data to 
guide in understanding how non-lincaritics can affect current 
interpretations of the 2dF power spectrum; A model fit to 
halo occupation number, based on the 2dF correlation func- 
tion, can be found in Magliocchetti & Porciani (2003) . In the 
middle and bottom panels of Fig. 4, we show the convolved 
redshift space power spectrum as a function of parameters 
in the halo description, mainly mass slope of the halo occu- 
pation number and the normalization of the virial equation 
that determines galaxy motions within halos. For compari- 
son, we also show the linear prediction, again convolved with 
the 2dF window function. One should note the lack of os- 
cillatory features in these power spectra when compared to 
those in Fig. 3. This is due to the averaging over the 2dF 
window function and leads to a reduction in amplitude of 
small scale features (e.g., Elgaroy, Gramann & Lahav 2002). 

As shown in Fig. 4, the redshift space galaxy power 
spectrum is such that the transition to non-linearity comes 
with a reduction of power, relative to the scaled linear power 
spectrum with a constant bias. This is opposite from the case 
associated with the real space power spectrum where the on 
set of non-lineaxity is always associated with an increase 
in power relative to the linear power spectrum scaled by 
a constant bias. Another way to state this observations is 
that the bias factor for galaxies in redshift space, relative 
to linear clustering, decreases from a constant value at large 
scales when non-linearities become important. As one moves 
to much small scales [k > 10 Mpc~^), the non-linear redshift 
space power spectrum dominates over the linear description 
and the bias factor increases back to a higher value (see. 
Fig. 3). This is consistent with calculations of the redshift 
space bias by Seljak (2001). As shown in Fig. 4, we found 
the reduction in non-linear power to be present when certain 
parameters related to the halo model description is varied 
and we consider this to be a general feature of the redshift 
space power spectrum. 

In general, our calculations indicate that the linear as- 
sumption is only valid out to ~ 0.1 Mpc^^ both in the case 
of the real and redshift space power spectra and proper mod- 
eling of these data beyond such a scale must include some 
accounting of non-linearities. If measurements beyond 0.1 
Mpc"'^ is used for cosmological purposes, with a constant 
bias factor, one can potentially misidontify the drop in power 
due to non-linearities as due to a cosmological effect, for ex- 
ample, due to a massive neutrino since massive neutrinos are 
expected to damp power at small scales. In the next Section, 
we address this issue in detail using the Fisher matrix based 
approach to quantify how non-linearities affect cosmologi- 
cal parameter determination and what improvements can 
be introduced to avoid potential biases. 




0.0 0.1 0.2 

k (Mpc') 



Figure 4. The redshift space power spectrum of the 2dFGRS 
(data points froiii Pcrcival ct al. 2001). The solid lines show the 
redshift space galaxy power spectrum convolved with the 2dF 
window function. Note the lack of baryon oscillations features 
due to smearing of modes over the finite width of the filter. For 
comparison, we also show the linear (dot-dashed line) and the 
redshift space dark matter (dotted lines) power spectra, again 
convolved with the 2dF window function. The middle and bottom 
panels are variations on the predicted redshift space galaxy power 
spectrum as a function of a, the mass-slope of the halo occupation 
number, and 7, the normalization of the virial equation. In the 
middle panel, we assume 7 = 1 while the bottom panel fixes a to 
be 0.8. In both case. No is taken to be 0.7. 
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Figure 3. The power spectrum of dark matter (dotted lines) and galaxies (solid lines) in the redshift space. For comparison, we also 
show the real-space power spectra which is lower than the redshift space at large scales and higher at non-linear scales. The increase in 
power at large scales is associated with bulk motions of halos while virial motions within halos lead to a suppression of power associated 
with the 1-halo term of the halo description. 



3 NON-LINEARITIES AND COSMOLOGY 

To quantify the extent to which cosmological parameters can 
be affected by assumptions in the literature and to calculate 
expected errors on cosmological parameters and associated 
biases, we make use of the Fisher information matrix: 



9^ InL 

dpidpj 



(17) 



whose inverse provides the optimistic covariance matrix for 
errors on the associated parameters (e.g., Tegmark et al. 
1997). In Eg. 1171 L is the likelihood of observing data set 
X, in our case the power spectrum, given the parameters 
pi . . .pn that describe these data. Following the Cramer- Rao 
inequality (Kendall & Stuart 1969), no unbiased method can 
measure the ith parameter with standard deviation less than 
(Fii)~^/^ if all other parameters are exactly known, and less 
than [(F~''")ii]^/^, from the inverse of the Fisher matrix, if 
other parameters are to be estimated from the data as well. 

For this discussion, we first focus on the binned real 
space galaxy power-spectrum and return to the case involv- 
ing the redshift space power spectrum as part of the dis- 
cussion. The binned power spectrum can be written as an 
integral average over a shell in k-space centered around some 
wave number ki with width Ak: 



(18) 



where Vi = Ank^Ak is the shell volume and V is the total 
survey volume. The full covariance related to measurement 
of this binned power spectrum is (Scoccimarro, Zaldarriaga 



& Hui 1999; Meiksin & White 1999; Cooray & Hu 2001) 



1 

V 



where 



d-^k f d-^k 



(19) 



(20) 



is the non-Gaussian contribution to the covariance of 
the power spectrum. This non-Gaussian contribution de- 
pends on parallelogram configurations of the trispectrum, 
T(ki, — ki,kj, — kj), and we make use of the 1-halo term of 
the halo approach (Cooray & Sheth 2002) to calculate it 
following Cooray & Hu (2001). This approach is a valid de- 
scription since in linear scales the covariance is dominated 
by the Gaussian variance while non-Gaussianities are only 
significant in non-linear scales where the 1-halo term domi- 
nates. We can write the relevant contribution as 



d^k 
'V~ 



dm n{m) 



d^k 

{iVgai...(jVgai -3)|m) 

"gal 



%al (fci I m) I ^ I Ugal (fcj j m 



where {Nga_i...{Ngai — 3)jm) is the fourth moment of the 
halo occupation distribution. We follow our previous de- 
scription related to the galaxy power spectrum and write 



(iVgai...iV(gal-3)|m) 



/3^(m)[2/3"(m) - l][3/3" 
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Figure 5. The ratio of modeled description of the galaxy power-spectrum to that of the total power spectrum as expected under the 
halo model. The case illustrated here is same as the one shown in Fig. 2. With the dotted line we show the linear power spectrum; as 
shown, the simple linear description with a constant bias lead to an underestimate of true power at the level of 20% when k ~ 0.2 to 
0.3 Mpc-i. This can be compensated with the addition of a power-law contribution. The ratio related to the linear plus a power-law is 
shown with a long-dashed line. The agreement is essentially exact out to 0.3 Mpc~^, but results in an overestimate of power at deeply 
non-linear scales of 1 Mpc~^ and at the level of 30%. For reference, with dot-dashed lines, we show the 1-halo contribution and the 
related power-law contribution that attempts to model the non-linear aspect related to the power spectrum. 



2] (A'gai)*, where departures from the Poisson model, de- 
noted with the function /3(m) in equation |3 is related to 
the adopted binomial distribution following Scoccimarro 
et al. (2001). In equation 1191 we allow for the shot-noise 
contribution to the power spectrum variance by redefin- 
ing P'i Pi + P^^ such that the power spectrum is av- 
eraged over the shell with the shot noise contribution of 
P^^ = l/fig, where fig = Ng/V is the number density of 
galaxies in the survey. In the Gaussian limit, the covariance 
only results in a diagonal error matrix with no correlations 
between band powers. This variance is given by the well 
known formula (Peacock & West 1992; Blake & Glazebrook 
2003) 



V 47rfc2Afc 



Pgal(fc) + - 



(22) 



Using the full covariance, the Fisher matrix for the 
galaxy power spectrum can be written as (Cooray & Hu 
2001) 



E 



dPa 



-idP^ 
dpi3 



(23) 



where, for example, pa denotes a cosmological parameter 
of interest. In addition to parameter errors. Fisher matrix 
approach also allows a useful technique to quantify the bi- 
ases associated with parameter estimates when the model 



description departs from the observations. For example, in 
the present case, the galaxy power spectrum is usually mod- 
eled with a combination of the linear power spectrum and a 
scale-free bias out to a certain wavenumber, beyond which 
no information is used for cosmological purposes. We can 
quantify the extent to which this is valid by studying pa- 
rameter biases such that if Pa = pa + 5pa where pa is the 
true value of the parameter a and 5pa is the resulting bias 
due to the wrong description of the power spectrum (in this 
case, linear power spectrum times constant bias), then 



dpp 



(24) 



where Pi and Pi{k) represent the true and erroneous binned 
galaxy power spectra at wave number ki, respectively (e.g., 
Huterer 2002). 

To calculate errors, we assume a galaxy redshift sur- 
vey consistent with the SDSS and take Vsdss = 10^ h""^ 
Mpc^ and Ng = 5 x 10^ such that the galaxy density is 
~ 5 X 10""^ h^ Mpc~^. For these assumed survey parame- 
ters, the non-Gaussian contribution to the covariance of the 
binned galaxy power spectrum is relatively minor and re- 
sults in a relative increase, when compared to the Gaussian 
covariance, by a factor of 1.4, at fe ~ 0.25 Mpc~^ and a fac- 
tor of 2 at fc ~ 0.5 Mpc~^. Note that if the number density 
of galaxies were to be higher, such that the shot-noise is re- 
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Figure 6. The la error ellipses for determination of Clmh^ and the neutrino mass. These errors follow from Tables 1 and 2. The solid line 
shows the constraint expected from the 3 year WM AP dataset (including polarization information) . The dotted lines show the expected 
improvement when galaxy clustering information is added with a cut off at A; < 0.1 Mpc~^ and using information out to 0.5 Mpc~^ 
with an accounting for non-lincaritics (sec text). The dashed lines show the same situation but with Planck as the source of information 
related to CMB. The galaxy clustering survey considered here is equivalent to that of the SDSS. 



duced at small scales, then the relative importance of the 
non-Gaussian contribution to full covariance increases. The 
power spectrum is assumed to be that modeled from the 
halo approach with a = 0.8 and A^o = 0.7. Wc rcnormal- 
ized the power spectrum such that the bias factor is unity 
at large scales. We then model the galaxy power spectrum 
through a combination of the linear power spectrum and the 
constant bias, again assuming to be unity. The parameter 
errors and associated biases are tabulated in Table 1 for a 
ACDM cosmological model with a dark energy equation of 
state w and allowing for massive neutrino contribution to 
the mass density of the universe. We tabulated errors and 
biases as a function of the maximum k value, fee, to which 
galaxy power spectrum analysis is considered to be linear. 
At typical large scales, with kc =0.1 Mpc~^, the biases are 
negligible and we recover parameter errors usually quoted 
in the literature by the Fisher matrix based analyzes sim- 
ilar to the one performed here (Eisenstein, Hu & Tegmark 
1998; Wang, Spergel & Strauss 1999). In Table 1, we also 
tabulate errors expected from the galaxy power spectrum 
alone and the galaxy power spectrum combined with final 
WMAP data (middle rows) and with Planck data (bottom 
rows). The CMB information also includes polarization data 
in addition to temperature. 

As one includes information out to smaller scales, or 
large fee values, from the galax:y power spectrum, certain pa- 
rameters, such as Qbh'^ and m^, become better determined 



than out to an early cutoflf, say at fee = 0.1 Mpc~^. This 
can be understood by the fact that the small scale regime of 
the galaxy power spectrum captures information related to 
these parameters through additional oscillatory features and 
the damping of power by massive neutrinos. This, however, 
comes at the expense of increasing biases in the parame- 
ter determination since in this regime, non-linearities enter 
the galaxy power spectrum as the 1-halo term becomes im- 
portant and starts to dominate when k > 0.6 Mpc^^ (see. 
Fig. 1). Thus, the simple description based on the linear 
power spectrum and the scale-free bias breaks down when 
kc > 0.1 Mpc^^ though there is significant amount of cos- 
mological information beyond this cut off scale. Note that in 
Table 1 we do not consider fee values above 0.5 since for the 
survey parameters we have considered, the shot-noise term 
Eissociated with the finite number density of galaxies starts 
to become important. For future surveys that attempt to 
measure an increased number of redshifts, it is clear that one 
gains further cosmological information by moving to smaller 
scales. 

While cosmological parameter estimation from data be- 
low ~ 0.1 Mpc~^ is unaffected by non-linearities, it is useful 
to consider possibilities that can return information at wave 
numbers above this value. This is due to the fact that cur- 
rent surveys, both 2dFGRS and SDSS as well as previous 
surveys such as IRAS PCSZ, probed to much smaller scales 
then implied by this arbitrary cut off. In the case of the real- 
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space galaxy power spectrum, as shown in Fig. 1 (bottom 
panel), we consider the possibility that the onset of non- 
linearity can be modeled by a power-law contribution when 
wave Numbers out to 1 Mpc~^ is considered. In fact, for a 

good approximation, the 1-halo term can be modeled as a 
power-law in the mildly non-linear regime with 

Pi^(fe) =^nifc"°' k<ki (25) 

where ki is the outer wave number to this approximation 
may be valid. This description includes three new parame- 
ters involving the normalization of the non-linear part of the 
power spectrum, Ani, the slope, ani, and the scale at which 
the power-law breaks, ki . Note that the power-law behavior 
remain the same when one considers different models of the 
halo occupation number (Fig. 2). 

The extent to which this approximation may be valid 
can be considered based on the halo model description of 
the non-linear galaxy power spectrum with the addition 
of a power-law contribution. While one can assume the 
halo model generated power spectrum as the expected non- 
linear power spectrum for galaxies, in parameter estima^ 
tion, one can consider the model fitting procedure based 
on the approximation that the non-linear power spectrum 
is P{k) = b'^(k)P^''^{k) + Anife""'. The last two parameters 
account for the non-linearity and, more importantly, the 
power law aspect of the 1-halo term. Since at large scales 
Wgai(fc|m) — > 1, one expects the normalization to be given 
by 

A.,= [dmnjm) (^^^.(iV.a. - l)|m) ^ ^^^^ 

J "gal 

and the power-law ani ^ (in Fig. 1, Qni ~ 3 since we 
plot the logarithmic power spectrum, or k''' P (k) / 2ii^ ) . The 
departure from this expected power-law captures the profile 
shape. Thus, if these two parameters can be extracted from 
the data as well, one not only obtains information related to 
cosmology, but some aspect related to how galaxies populate 
dark matter halos. The numerical value of A^\, for example, 
can aid in constraining model independent measurements of 
the halo occupation number (Cooray 2002b). 

The extent to which this approximation may bo valid 
is summarized in Fig. 5, for the same case shown in Fig. 2. 
Here, wc plot the ratio of the model power spectrum to that 
of the non-linear case as calculated by the halo model. As 
already discussed, the description that P{k) — 6'^(fc)P''"(fc) 
breaks down rapidly such that at a wavcnunibcr of 0.1 
Mpc^^ one underestimates the true non- linear power by a 
few percent while at 0.3 Mpc~^, this underestimate is at 
the level of 25%. Instead of P{k) = b'^{k)P^'"'{k), when 
one includes the additional power-law such that P{k) = 
h^{k)P^"^{k) + Anifc""', the difference is below a few per- 
cent out to a wavenumber of 0.3 Mpc"'^ suggesting that 
the transition is well modeled with the single power-law ap- 
proximation. At a wavenumber of 0.5 Mpc"'^, the power- 
law assumption leads to an overestimate of power at a level 
of 5% and this overestimate is at the level of 30% when 
fe ~ 1 Mpc"'^, deep in the non-linear regime. While the sin- 
gle power-law assumption clearly breaks down at these small 



scales, the improvement over the simple linear assumption 
is significant when k ~ 0.3 Mpc~^. 

While we have only considered a specific case with cer- 
tain parameters for the halo occupation number in Fig. 5, 
we expect the situation to be the same with different de- 
scriptions. The extent to which one departs, however, is 
model dependent as wc have shown in Fig. 2, especially 
with respect to the slope of the halo occupation number- 
halo mass relation (Fig. 2 top panel). The differences move 
the non-linear scale from 0.05 to 0.2 Mpc^^ suggesting that 
the importance of the power-law varies over the same scales 
for different samples of galaxies. The extent to which the 
power-law approximation works will similarly vary and we 
note that depending on a, this cam be between 0.3 and 0.5 
Mpc"'^ while departures are at the level of 10% to 50% at 
1 Mpc^'^; if the power-law approximation is not included, 
the departures are substantial (> 100%). Regardless of the 
exact variations, such as on the extent to which the power- 
law assumption overestimate power when fc ~ 1 Mpc~^, the 
power-law assumption provides a useful addition to the sim- 
ple linear power spectrum description, with a constant bias, 
in mildly non-linear scales of order 0.2 to 0.4 Mpc~^. 

Note that the power-law assumption is independent of 
selection effects and biases related to observational aspects. 
This can be understood simply as following. The observed 
total power spectrum is simply a convolution of the real 
power spectrum time a window function that captures ob- 
servational issues (Eq. 16). Since the true power spectrum 
can be described better with a combination of the linear 
power spectrum and a power-law correction, in the mildly 
non-linear scales, the observed power spectrum will remain 
essentially described by the same combination with the fil- 
tering aspects related to observations applied. Thus, regard- 
less of how observations are done or the galaxy power spec- 
trum is measured, we expect our description to be accurate 
to the level we can also describe the true underlying power 
spectrum with the same assumptions. 

We can further understand the extent to which this 
model is valid by considering its ability to extract cosmologi- 
cal parameters. In Table 2, we list expected errors on cosmo- 
logical parameters in addition to these two new astrophys- 
ically oriented parameters . While cosmological parameter 
estimates either remain same as in Table 1 or worsen slightly, 
we note there is a sharp reduction in biases associated with 
parameter measurements when power spectrum information 
out to wave rmmbcrs of 0.5 Mpc^^ is included. This is sig- 
nificant since the assumed power-law behavior around the 
onset of non-linearity allows one to measure cosmology from 
previously ignored non-linear scales such that no biases are 
introduced. In Fig. 6, as an example, we show the measure- 
ment errors for two parameters {Q.mh? and mi,) by marginal- 
izing over all other parameters. We show expected errors in 
combination with WMAP and Planck and the improvement 
when information out to 0.5 Mpc~^ is included with two ad- 
ditional parameters to be determined from data. Including 
information out to 0.5 Mpc~^ improves parameter errors, for 
example in the case of neutrino mass, by roughly a factor of 
2. 
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So far, we have only considered the real space galaxy 
power spectrum. A similar situation exists for the redshift 
space power spectrum, though, as discussed in Section 2.2, 
the onset of its non-linearity is associated with a decrease 
in power, relative to linear power spectrum times a constant 
bias, instead of an increase. To model this behavior, one can 
introduce an additional damping term associated with non- 
linear pairwise velocity to the usual formula involving the 
linear power spectrum (Ballinger, Heavens & Taylor 1996; 
Hatton & Cole 1998; Kang et al. 2002). This leads to a 
modification to equation |n| 



[1 + 



[l + fcV^cr2/2]'''' 



(27) 



such that the angular averaged redshift space power spec- 
trum takes the form of 



6(fc2a2/2)^ 



+ 



(/3 - fcVV2) arctan(^fc2cr2/2) 



(28) 



where /3 = f{Qm)/bg, as before, and bg is the scale free 
large scale galaxy bias. When ct — > 0, this formula is equiva- 
lent to the standard formula used in cosmological interpre- 
tations with the linear power spectrum scaled by the fac- 
tor (1-1- 2/3/3 -|- 1/5/3^ )6g. Note that such a correction only 
accounts for the increase in power at large scale associated 
with bulk motions. At small scales, with wavenumbers ~ 0.1 
Mpc~^ and larger, it is necessary to account for suppression 
of power related to virial motions with dispersion a, and we 
suggest the use of equation 1281 In order to account for the 
1-halo part of the non-linear redshift space power spectrum, 
especially for an analysis of data beyond a wavenumber of 
~ 1 Mpc~^, one can introduce again a power-law, following 
the suggested approach for the non-linear real space power 
spectrum. At mildly non-linear scales, the approach sug- 
gested here leads to an additional parameter, a, which can 
be simultaneously determined from data in additional to 
cosmological estimates. Currently, the well studied redshift 
space power spectrum comes from the 2dF survey and we 
plan to return to these issues as part of a combined mod- 
eling effort on non-linearities in this power spectrum based 
on cosmology as well as simplified corrections based on the 
halo approach. 



4 SUMMARY 

The galaxy power spectrum is now a well known tool of pre- 
cision cosmology. Ongoing surveys such as 2dF and SDSS 
have already provided a wealth of information related to cos- 
mological parameters based on clustering of galaxies in the 
presumed linear regime, though measurements span down 
to much smaller scales than used for cosmological purposes. 
We have shown here that there is considerable amount on 
information in the regime where clustering transits from lin- 
ear to non-linear. The extraction of this information, how- 
ever, cannot be done alone with the usual approach involv- 
ing the linear power spectrum and a constant scale-free bias. 



We have investigated the non-linear behavior of the galaxy 
power spectrum in both real and redshift space based on 
the halo approach and have suggested that the on set of 
non-linearities can be modeled as a simple power-law. The 
normalization and the power-law slope provide certain in- 
formation related to halo occupation description in the halo 
models of galaxy clustering. We have shown that significant 
improvements can be made when non-linear effects are taken 
into account through this power-law contribution with two 
additional parameters to be determined from the data si- 
multaneously with cosmological parameters. The suggested 
approach can easily be implemented in current cosmological 
studies of galaxy clustering and to extract most information 
from galaxy clustering data to date. 
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